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ABSTRACT 

We present results based on Chandra observations of a large sample of 129 hot galaxy 
clusters. We measure the concentration parameter C200, the dark mass M200 and the 
baryonic mass content in all the objects of our sample, providing the largest dataset of 
mass parameters for galaxy clusters in the redshift range z = 0.01 - 1.4. We confirm 
that a tight correlation between C200 and M 2 oo, c oc M£ ir / (l + z) b with a = -0.56 ± 0.15 
and b =0.80± 0.25 (68 per cent confidence limits), is present, in good agreement with 
the predictions from numerical simulations and previous observations. Fitting the mass 
profile with a generalized NFW model, we got the inner slope a, with a — 0.94 ±0.13. 
Finally, we show that the inner slope of the density profile, a correlates with the 
baryonic mass content, namely a is decreasing with increasing baryonic mass 
content. 

Key words: cosmology-theory-large scale structure of Universe-galaxies-formation 



1 INTRODUCTION 



The problem of dark matter halos formation has 
a long histo ry goin g bac k to th e sem inal paper of 
iGunn fc Gottl (|l972r ) and iGunnl l| 19771 ). who studied 
the density profile formation using the collapse of a 
spherical perturbation in an expanding background. The 
quoted papers and several following analytical models 
(e.g.. iFillmore fc Goldreichl 1 1984 iBertschingel 1 19851 ; 
iHoffman fc Shahamlll985l ), found that density profiles are 



Bullock et ail l200ll IPower et all 12001 and iNavarro et~ai] 
2004 INavarro et al I I2OI0F found that the spherically aver- 
aged density profiles of the N-body DM halos are similar, 
regardless of the mass of the halo or the cosmological 
model. The NFW profile is given by: 



p{r) = 



f>(> 



(1) 



described by power-laws in all the ra dius range. More 



1992!, Subramanian et al. 200d; Del Popolo et al. 




2000; 


El-Zant et all 2001 . El-Zant et al. 2004 Hiotelis 


2002; 


Le Delliou & Henriksenl 20031; 


Ascasibar et al. 




2004; 


Williams et all |2004|; Tonini et al 


2009; Ascasibar et al. 


2007f) showed that the profile is not a power-law, similarly 



r/r s (l + r/r s ) 2 r/r s (l + r/r s ) 2 

where p cr itic is the critical density of the Universe at the 
cluster's redshift 2, and A c is the virial overdensity. The 
scale radius r s is connected to the virial radius r V i r through 
the concentration parameter c, c — r V i r /r^\. At small scales, 
the logarithmic density slope, that is also known as inner 
slope, is given by 

dlogp 



Navarro et"aH i 19971)); Navarro-Fre n k- White (NFW)) , 
19981 ; IJing fc Sutd l2000l : iKlypin et all l200ll 



dlogr 

The quoted profile diverges as p oc r _1 
part, and at large radii behaves as p oc r~ 



(2) 

in the inner 
3 . The inner 



Moore et al 
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Since r v i r is difficult to determine observationally, its value 
is often approximated by the radius within which the average 
density is greater than the critical density by a specified factor 
(e.g., 200). In the following of the paper, r v i r is identified with 
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slope is even steeper in I Moore et al.l il 19981) profile (p oc 



-1.5\ 



More recent 
7\U. 



Power et al] 120031. 



si mulations (e.g., ^ 
Havashi et al.ll2004l andlNavarro et al.ll2004l . iNavarro et aD 
20ld ; IStadel et all l2009h showed that density profiles are 



better fitted by the Einasto profile, which becomes shallower 
towards the centre of the halo. 

Unfortunately, N-body simulations predictions are 
not supported by observations. Observations of the in- 
ner part of density profiles of dwarfs galaxies and 
LSBs are characterized by a core-like structure (e.g. 



Flores & Primackl Il994; 


Moore! Il994l; Ide Blok & Bosma 


2002 


;lde Blok et alj 2003; 


Gentile et al. 2004, Salucci et al. 


2007 


; Kuzio de Narav et al.1 20081. Kuzio de Narav et al. 


2009 


), and a similar problem is evidenced when studying 



the clusters of galaxies inner profile. Density profiles of clus- 
ters has been studied through X-ray observations, strong 
and weak lensing. X-ray temperature measurements give 
information on clus ter structure in the range 500-50 kpc 
(Brad aFet al.l [2008 ). At smaller radii, temperature deter- 
mina tion is limited by instru mental resolution or substruc- 
ture |Schmidt fc Alien! [20071 . hereafter SA07). X-ray mea- 
surements are also limited by "cooling flows" presence, and 
the breaking of assum ption of hydrostatic equilibrium (see 
lArabadiis etail 120041 '). The inner slope calculation through 
the use of X-ray observations broug ht to discrepant v alues 
( e.g.. 0.6 (lEttoriet al.ll2002f ). 1.2 (|Lewis et al.ll2003l ), 1.9 
jArabadiis et al.ll2002T v 

Another technique used to study the DM distribution in 
clusters is gravitational lensing. Weak lensing of background 
galaxies is used to reco nstruct the m ass distribution in the 
outer parts of clusters (|Mellierl 1 19991 ). The resolution that 
can be achieved is able to constrain profiles inside 100 kpc. 
Strong lensing is used to study the DM distribution in the 
inner parts of clusters, it has a typical sensitivity to the 
projected mass distribu tion inside ~ 100 — 200 kpc, with 
limits at ~ 10-20 kpc i|Gavazzill2005l ; iLimousin et alJl200St 
Alexandrov et al. 2011; Tsvetkova et al. 2009). 

Discrepant results ha ve been sometime o btained when 
using the lensing method. ISmith et al.l l|200ll ). found a > 1 
at 1% of r vir by studying the tangential and ra dial arcs of 
A383. A m uch smaller value a w as obtained bv lSand et al.l 
l|2004l) andU cwman et all l|201ll) , for the same cluster, us- 
ing lensing an d through the aid of stellar kinematics of the 
central region. iTvson et all il 19981) found a = 0.57 ± 0.024 
for CI 0024+1654, while iKneib et al.1 (|2003l ) found that a 
NF W profile fits the p rofile in the radius range 0.1-several 
r,„>. ISand et al.l l|2002h f ound a cored p r ofile w ith a = 0.35 
for M S2137.3-2353, while iGavazzi et al.l (|2003l ) and lGavazzil 
|2005l) concluded that the precise value of the slope depends 
on the mass-to-light ratio of the Brightest Cluster Galaxy 
(BCG). 

In summary, N-body simulations are not always in 
agreement with the inner slopes of dwarf galaxies, LSBs 
and clusters of galaxies, and in the case of clusters of 
galaxies , observations may even disagree for the very same 
clust er jGa vazzi 200 51 ; ISmith et al . 2005; Z appacosta et all 
20061; SA07; iBradac et al.l 12008? ILimousin et al.l l2008t 



Umetsu fc Broadhurstl l2008). 



Several could be the reasons of the discrepancy (e.g., 
(a) different definition of the slope, which sometimes refers 
to the DM and sometimes to the total mass; (b) use of ob- 
servational techniques with different/limited dynamic range 



in radius; (c) not taking into account the stellar mass of 
the BCG). In order to obtain more sure constraints on 
the central part o f the density profiles, one should use 
combined methods jMiralda-Escud e 1995; K neib et al.ll2003l 
(weak+s trong lensing) ; Bradac et al. 20051 (weak+strong 



lensing) ; lMahdavi et al.ll2007l (X-rav+weak lensing)), or bet- 
ter constraints on the central part of the density profiles can 
be obtained through stellar kinematics of the central galaxy 

(~ 1 — 200 kpc region). 

In a series of paper s, ISand et al.1 (|2002l ). ISand et al.l 

l|2004l ). ISand et al.1 (|2008l ). studied the clusters MS 2137-23; 
A383; A963; RXJ1133; MACS 1206; A1201, separating the 
contribution to the halo coming from the DM from that 
coming from the baryonic and stellar mass of the BCG. 
They found a p r ofile fl atter than a < 1 except for RXJ1133. 
INewman et al.l (1201 ll) found a < 1 (95 per cent CL) for 
A383. iNewman et al.1 (|2009l ) presented a detailed analysis 
of DM and baryonic distributions in A611, finding a slope 
a < 0.3 (68 per cent CL). 

Summarizing, at least some clusters of galaxies have 
inner density-profile slopes shallower that those obtained 
in N-body simulations, in agreement with what happens, 



galaxies and LSBs (e.g. Flores & Primackl 1994 


; Moordll994l; 


de Blok & Bosmal |2002|; Ide Blok et all 120031; 


Gentile et al.l 


20041. Salucci et all 12007 


; Kuzio de Narav et al.l 20081. 


Kuzio de Narav et al.l I2OOS 


). We want to recall here that, 



similarly to Sand's result concerning the quoted clusters, 
the galaxies NGC 2976, NGC 4605, NGC 5949, NGC5963 
and NGC 6689, have inner slop es going from very flat to 
cuspy, and Ide Blok et al.l ()2008h , using a sample from The 
HI Nearby Galaxy Survey (THINGS) found that the best 
fit to rotation curves, and then the inner slope a of their 
density profile, depends on their mass. 

Several papers have shown the fundamental role of 
baryons in shaping the density profiles of structures. Differ- 
ent processes have been pointed out capable of flattening 
the inner density profile, transferring energy from stellar 
baryons to the dark matter, heating it and lowering the 
central dark matter density dMilosavlievic fc Merrittll200ll; 



El-Zant et al.l l200ll. |E1-Zant et al.1 l200ll; IWeinberg fc Katz 



Loeb fc Peebles! l2003l; lGao fc White! 120071; iMa et al 



iMcMillan fc Dehnenl p005l; iTonini et al.l 
-'- ' 1 2006; iRomano-Dfaz et al 



2006 



2008, 



2002 
2009 

Mashchcnko et al 
[ Romano- Diaz et al.ll2009l ; lDel Popo lo 2009; Governato et al.l 
120101 ; Kulinich et al. 2012). 

In the present paper, we use Chandra X-ray data to 
study, similarly to SA07, the properties of 129 dynamically 
relaxed galaxy clusters. We selected the sample with red- 
shift range from 0.01 to 1.4. with the aim to recover their 
total and gas mass profiles and analysis of the measured 
distribution of C200, M200 and baryonic mass content. We 
determine the baryons content of each cluster subtracting 
the DM from the total mass. Even if the results concerning 
the density profile of clusters are similar to those of SA07, we 
point out that fitting the density profile with a generalized 
NFW model 



p(r) 



Po 



(r/r s ) a (l+r/r a ) 3 - a 



(3) 



the slope a is correlated with the baryonic mass con- 
tent of the cluster, in agreement with several studies 
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Figure 1. The distribution of 
sphere. 
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Figure 2. The distribution of galaxy clusters in our sample. 



jRicottil l2003l: ISubramanian et al.l I2000I; ISimon et al.ll2005l; 
Cen et all I2004I: IRicotti fc Wilkinson! l2004l: IWi liams et all 
2004 IRicotti et all I2007I ; iDel Popold I2009I . iDel Popolol 



2012J). In the paper the assumed value of Hubble constant 
is Ho = 73fcms _1 Mpc _1 and S7 m = 0.27. The outline of our 
work is the following. In Section 2, we describe the X-ray 
galaxy clusters sample of Chandra observations compiled 
by us to recover the baryonic and total mass profiles with 
techniques presented in Section 3. In Section 4, we present 
a preliminary discussion of the main results. We investigate 
the C200 — M200 and a — Mt relations. We summarize our re- 
sults and draw the conclusion of the present study in Section 
5 and 6. 



2 SAMPLE AND DATA PROCESSING 

The sample has been built to cover a wide range in red- 
shift from 0.01 to 1.4, constituted by 129 galaxy clusters. 
All clusters in the sample have a regular X-ray morphology, 
indicative of a relaxed state and allowing reliable determina- 
tion of the total mass profile through the hydrostatic equi- 
librium equation. The observations results are presented in 
Tab. [T] The preliminary processing of the Chandra data, 
which where obtained using the Advanced CCD Imaging 
Spectrometer (ACIS), included: the Chandra pipeline pro- 
cessing of the list of events; the re-processing with the CI AO 
software package; improving the rejection of cosmic ray 
events by VFAINT mode; cleaning the data in format of 
recommendation by the Chandra X-ray Centre (to remove 
periods of anomalously high background). 

The Tab. 1 consists of the name of clusters, the red- 
shift, the number of observation, the exposition time (be- 
fore cleaning) , instrume nt, and the values of c olumn density, 
which were described bv lDickev fc Lockmanl ( |l990h . and co- 
ordinates (RA, DEC) (which were taken from NEE0). It is 
important to note that our list of clusters includes a few 
clusters with merger's properties, for example, Abell 2744, 
but in such case we have used only one component on the 
image for data reduction. 

As we see in Fig[2] the most of clusters has the redshift 
in the range 0.01 - 0.3. We added to our sample MACS and 
RCSJ clusters as ones of the most distant objects which are 
observed by Chandra. The distribution of the studied galaxy 



http:/ /nedwww. ipac.caltech.edu 



clusters at the celestial sphere is given in Fig. [T] We used in 
our re search the sample of the clusters used by en et al.l 
(2002) to study the evolution of the X-ray gas mass fraction 
and constrain cosmological parameters. To minimize system- 
atic scatter and to allow the most precise test of the CDM 
model predictions, we used only highly relaxed clusters. 

The data reduction was performed with the Chandra 
CIAO v.4.2 including CalDB 4.5.1 for maps and calibra- 
tion. The main data pro cessing steps w ere provided using 
techni q ues discussed by Babvkl (120121') . iBabvk fc Vaviloval 
(|2012l ). lBabvk et all (|2012l ). IBabvk et all (|2012h in our pre- 
vious research. The quoted steps are the following: a) re- 
moving the point sources; b) determining the X-ray peaks; 
c) centering the X-ray peacks with different width to provide 
the same number of X-ray photons in each annulus; d) ex- 
tracting the annular spectra. The average of the outermost 
radius is about 1 Mpc. We have generated ARF and RMF 
files using the mkarf and mkrmf command in CIAO tool. 

We used Xspec (version 12.6)(Arnaud 1996) to analyze 
spectra, the MEKAL0 code (Kaastra & Mewe 1993) and 
WABqj for fitting the data and Fe-L calculations (Liedahl 
et al. 1995). The content of heavy elements in this model was 
taken close to the solar level (Z=0.3) and was frozen. The 
energy range 0.5-7.0 keV was used during fitting, and the 
spectra were grouped to have at least 50 counts per spectral 
bin. 

Background spectra were extracted from the blank-field 
data sets which is available from the Chandra X-ray centre, 
and this technique was used only for nearer galaxy clusters 
(0< z <0.3). For more distant clusters, the background spec- 
tra were extracted from the region which has the same size 
like in region for source. Usually, we used observations which 
were made on ACIS chips 0,1,2,3 and 7, because these chips 
are most accurate in calibration, although ACIS 4,5 and 6 
were also used. All background spectra were cleaned as done 
with the main previous spectra. 

We have used DSDEPROJ rou tine method 
i Sanders fc FabiarJ[2007l . iRussell et al.| [2008) to process the 
annular spectra in order to determine X-ray temperature 
and other parameteres. After determining the temperature 
of clusters we built the surface brightness profile for each 



3 MEKAL is the model which describes an emission from hot 
diffuse plasma (ICM) 

4 Th e WABS is a parameter w hich describe the galactic absorp- 
tion jPickev fc Lockmanlli99rj) 
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Table 1. Summary of the Chandra observations. Column list the target name, redshift, number of ID observation, net exposure 
time, detector used, column density and coordinates (taken from NED) from the X-ray centres used in the analysis. 



Cluster z ObsID t<, xp , Instrument n H , RA(J2000) DEC(J2000) 
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00 
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-40 03 00 
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5.58 


06 


26 


15 


.4 


-53 40 52 


A3395 


0.050 


4944 


22.17 


ACIS-I 


6.13 


06 


27 


31 


.1 


-54 23 58 


A3526 


0.011 


504 


32.12 


ACIS-S 


8.07 


12 


48 


51. 


,8 


-41 18 21 


A3558 


0.048 


1646 


14.61 


ACIS-S 


3.89 


13 


27 


54. 


,8 


-31 29 32 


A3562 


0.049 


4167 


19.54 


ACIS-I 


3.83 


13 


33 


31 


,8 


-31 40 23 
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Table 1. Continued. 



Cluster z ObsID texp, Instrument n H , RA(J2000) DEC(2000) 

ks 10 20 cm" 2 



A3571 


0. 


.039 


4203 


34.44 


ACIS-S 


3.91 


13 


47 


28.9 


-32 51 57 


A3667 





.055 


889 


50.96 


ACIS-I 


4.95 


20 


12 


30.1 


-56 49 00 


A3827 


0. 


.098 


7920 


46.17 


ACIS-S 


2.12 


22 


01 


49.1 


-59 57 15 


A3921 


0. 


.093 


4973 


29.76 


ACIS-I 


2.82 


22 


49 


38.6 


-64 23 15 


A4038 





.030 


4992 


33.97 


ACIS-I 


1.54 


23 


47 


37.0 


-28 07 42 


A4059 


0. 


.047 


897 


41.19 


ACIS-S 


1.10 


23 


56 


40.7 


-34 40 18 


AWM4 


0. 


.032 


9423 


75.46 


ACIS-S 


5.03 


16 


04 


57.0 


+23 55 14 


CLJ1226.9+3332 





.890 


3180 


32.12 


ACIS-I 


1.37 


12 


26 


58.0 


+33 32 54 


Coma (A1656) 


0. 


.023 


10672 


28.91 


ACIS-S 


0.93 


12 


59 


48.7 


+27 58 50 


IIIZw54 


0. 


0311 


4182 


23.76 


ACIS-I 


15.50 


03 


41 


17.6 


+15 23 44 


ISCS J1438+3414 


1. 


.41 


10461 


150.00 


ACIS-S 


2.01 


14 


33 


22.4 


+22 18 35.4 


MACSJ0011. 7-1523 


0. 


.379 


3261 


22.00 


ACIS-I 


2.05 


00 


11 


42.9 


-15 23 22 


MACSJ0159. 8-0849 





.405 


3265 


18.00 


ACIS-I 


2.08 


01 


59 


48.0 


-08 49 00 


MACSJ0242. 6-2132 


0. 


.314 


3266 


12.00 


ACIS-I 


2.69 


02 


42 


36.0 


-21 32 00 


MACSJ0429. 6-0253 





.399 


3271 


24.00 


ACIS-I 


5.25 


04 


29 


36.0 


-02 53 00 


MACSJ0647.7+7015 


0. 


.59 


3196 


18.85 


ACIS-I 


5.63 


06 


47 


45.9 


+70 15 03 


MACSJ0744.8+3927 


0. 


.697 


536 


19.69 


ACIS-I 


5.68 


07 


14 


51.8 


+39 27 33 


MACSJ1115. 8+0129 


0. 


.355 


3275 


16.00 


ACIS-I 


4.39 


11 


15 


53.3 


+01 29 47 


MACSJ1311. 0-0311 


0. 


.49 


3258 


15.00 


ACIS-I 


1.88 


13 


11 


00.0 


-03 11 00 


MACSJ1423. 8+2404 





.545 


4195 


113.40 


ACIS-I 


2.38 


14 


23 


48.3 


+24 04 47 


MACSJ1427.6-2521 





.318 


3279 


17.00 


ACIS-I 


6.18 


14 


27 


39.4 


-25 21 02 


MACSJ1720. 3+3536 


0. 


.391 


3280 


21.00 


ACIS-I 


3.37 


17 


20 


15.5 


+35 36 21 


MACSJ1931. 8-2635 


0. 


.352 


3282 


14.00 


ACIS-I 


9.11 


19 


31 


48.0 


-26 35 00 


MACSJ2129. 4-0741 


0. 


.589 


3199 


17.69 


ACIS-I 


4.84 


21 


29 


26.0 


-07 41 28 


MACSJ2229. 8-2756 





.324 


3286 


17.00 


ACIS-I 


1.34 


22 


29 


48.0 


-27 56 00 


MKW3s 





.045 


900 


58.03 


ACIS-I 


3.04 


15 


21 


51.9 


+07 42 31 


MKW4 


0. 


.020 


3234 


30.36 


ACIS-S 


1.90 


12 


03 


57.7 


+01 53 18 


MKW8 





.027 


4942 


23.45 


ACIS-I 


2.80 


14 


10 


43.1 


+03 27 11 


PKS0745-19 


0. 


.103 


508 


28.33 


ACIS-S 


42.70 


07 47 31.32 


-19 17 40.0 


RXCJ0043. 4-2037 





.292 


9409 


20.18 


ACIS-I 


1.54 


00 


43 


23.1 


-20 37 35 


RXCJ0232. 2-4420 





.283 


4993 


23.71 


ACIS-I 


2.42 


02 


32 


18.7 


-44 20 41 


RXCJ0307.0-2840 


0. 


.253 


9414 


19.16 


ACIS-I 


1.36 


03 


07 


01.1 


-28 40 30 


RXCJ0516. 7-5430 


0. 


.294 


9331 


9.64 


ACIS-I 


6.80 


05 


16 


35.2 


-54 16 37 


RXCJ0547.6-3152 


0. 


.148 


9419 


20.04 


ACIS-I 


2.02 


05 


17 


34.2 


-31 53 01 


RXCJ0605. 8-3518 





.141 


12899 


5.06 


ACIS-I 


4.47 


06 


05 


52.4 


-35 18 22 


RXCJ1131. 9-1955 





.307 


3276 


14.10 


ACIS-I 


4.46 


11 


32 


00.7 


-19 53 34 


RXCJ2014.8-2430 


0. 


.161 


11757 


20.18 


ACIS-S 


7.77 


20 


11 


49.7 


-24 30 30 


RXCJ2129. 6+0005 





.235 


552 


10.09 


ACIS-I 


4.29 


21 


29 


37.9 


+00 05 39 


RXCJ2337.6+0016 


0. 


.273 


3248 


9.31 


ACIS-I 


3.82 


23 


37 


39.7 


+00 17 37 


ZwCL1215 





.075 


4184 


12.22 


ACIS-I 


1.74 


12 17 40.60 


+03 39 45.0 


RCSJ0224-0002 





.778 


4987 


90.15 


ACIS-S 


2.92 


02 24 00.0 


-00 02 00 


RCSJ0439-2904 


0. 


.951 


3577 


77.17 


ACIS-S 


2.63 


04 39 38.0 


-29 04 55 


RCSJ 1107-0523 





.735 


5825 


50.12 


ACIS-S 


4.24 


11 07 22.80 


-05 23 49.0 


RCSJ1419. 2+5326 


0. 


61 


3240 


10.03 


ACIS-S 


1.18 


14 19 12.0 


+53 26 00 


RCSJ1620+2929 


0. 


.87 


3241 


37.13 


ACIS-S 


2.72 


16 20 09.40 


+29 29 26.0 


RCSJ2156+0123 





.335 


674 


45.76 


ACIS-S 


2.56 


15 47 34.2 


+26 38 29.0 


RCSJ2318+0034 


0. 


.78 


4938 


51.12 


ACIS-S 


4.13 


23 18 30.67 


+00 34 03.0 


RCSJ2319+0038 


0. 


.904 


5750 


21.18 


ACIS-S 


4.16 


23 19 53.00 


+00 38 00.0 


RXJ0439.0+0520 


0. 


.208 


527 


9.71 


ACIS-I 


10.50 


04 39 02.2 


+05 20 43 


RXJ0848.7+4456 


0. 


.570 


927 


126.74 


ACIS-I 


2.66 


08 48 47.2 


+44 56 17 


RXJ0849+4452 


1. 


26 


945 


128.45 


ACIS-I 


2.50 


08 53 43.6 


+35 45 53.8 


RXJ0910+5422 


1. 


.106 


2227 


84.20 


ACIS-I 


2.35 


09 10 45.36 


+54 22 07.3 


RXJ1113.1-2615 





.730 


915 


105.95 


ACIS-I 


5.47 


11 


13 05.2 


-26 15 26 


RXJ1221.4+4918 


0. 


.700 


1662 


80.13 


ACIS-I 


1.45 


12 21 24.5 


+49 18 13 



cluster. DSDEPROJ is the deprojection routine which 
assumes only spherical geometry, and solves some of the 
issues inherent to model-dependent deprojection routines. 
DSDEPROJ produces a set of "deprojected spectra" which 
can then be fitted by a suitable spectral model in Xspec. 
We have shown that this method does not generate the 
oscillating temperature profiles for multi-temperature 
clusters and produces a stable solution for an elongated 
cluster and clusters with breaks in temperature or densityQ 
For each ring we determined the temperature, kT, and 
the parameter norm ~ J n^nndV which is proportional to 



the electron, n e and hydrogen, tih, concentrations. Note, 
that the other parameters of model were fixed. Using these 
parameters, we have built the surface brightness profiles 
which fitted our theoretical surface profiles and then it was 
used to obtain the mass of each cluster (see below). The 
surface brightness profiles are measured in the 0.5-7 keV en- 
ergy band, which provides an optimal ratio of the cluster 
and background flux in Chandra data. 



5 The DSDEPROJ source code is available at www- 
xray.ast.cam.ac.uk/papers/dsdeproj. 



Note, that for lower redshift clusters in our sample, the 
statistical accuracy of the surface brightness at large radii is 
limited mostly by the Chandra field of view. 



6 Iu. Babyk, I. Vavilova, A. Del Popolo 



3 METHODS 

3.1 Galaxy clusters mass profile 

In order to calculate the total mass profiles of galaxy clus- 
ters, spherical symmetry and hydrostatic equilibrium of clus- 
ters was assumed. Similarly to SA07, we modeled the total 
mass profile as composed of three components: the dark mat- 
ter halo, the gaseous component and the optically luminous 
mass of the central cD galaxy. 

For the density profile of DM, we used the NFW model 
having a DM matter content within a fixed radius given by 



Af(< r) = 4ivp r s 



ln(l + —)- 



(4) 



The gravitational potential 
lated by the formula 

d4 _ C M{< r) 
dr r 2 



of dark matter can be calcu- 



(5) 



We made all the computations assuming the hydrostatic 
equilibrium condition of the X-ray emitting gas and wrote 
the hydrostatic equilibrium condition as follows 



V P = -Pa V ' 



(6) 



where P and p g are the gas pressure and the density, re- 
spectively. Since the density and pressure of hot gas have 
very low values, we can use the ideal gas law, P = P{1 , 
where fj, is molecular weight, and m p the proton mass. We 
then obtained the equation for unknown gas density distri- 
bution 



VPs 
Pg 



V' 



pm p 
kT g (r) 



(7) 



For the construction of a hot gas density scalar field of galaxy 
clusters, we had to calculate numerically the system of the 
following differential equations 

1 dp g pm p dcj){r) 

p g dxi kT g (r) dxi 

for gas distribution, where Xi = (a;, y, z) are the Cartesian 
coordinates, and then for total mass 



(8) 



M(< r) 



(9) 



In case of a single spherically-symmetric galaxy clus- 
ters, it is possible to calculate gas density a nd surface bright- 
ness p rofile analytically, like i n the w orks of lHumphrev et aD 
l|200rj) and lO' Sullivan et all (|2007h . for example. We per- 
formed the numerical integrations of Eq. ([§} and Eq. © by 
means of the Euler method. The integration step was cho- 
sen in order to obtain the maximum error of resulting X-ray 
brightness not bigger than 0.1 per cent. 



the model, we can obtain the gas density with precision of 
some constant factor a which depends on the integration 
boundary conditions of Eq. Therefore the real density 
p g cal of gas is equal to p r g eal = a-p s g lm , where p s g lm is the den- 
sity of gas obtained solving the previous differential equa- 
tions, that we indicate with sim standing for simulations. 
So, we can write p s g m = p g in Eq. ((8j). Thus the value of 
emission measure EM S i m from simulations can be found as 



EAIsim — 



j n e n p dV = 0.64/m 2 J (p S g lm ) 2 dV. (10) 



Using equation fllUp . we found the next expression 

norm alm = EM sim ■ 10" 14 /[4tt(I>a(1 + z))% (11) 

where Da is the angular distance to cluster. From the 
observational data we got the normalization parameter of 
the best-fit MEKAL model tiotytim ekal which was fit- 
ted by the norm a i m values from the simulations. Since 
norrriMEKAL parameter can be expressed by 



norrriMEKAL 



0.64 f (p r g eal ) 2 dV ■ 10- 



(12) 



from Eq. (j 101 Hip and Eq. (|12[1 , we have found the normal- 
ization factor a as 



norrriMEKAL 



(13) 



Therefore we able to obtain the field of hot gas density 
Pg eal — a ■ p g ' m and emission measure EM Tea i putting into 
\Mp7 al instead of p 3 g lm . 

For each cluster we made the series of numerical simula- 
tions of the surface brightness profile reconstruction. In our 
model we used two free parameters po and r s from NFW 
model. Sorting them out in the acceptable ranges, we se- 
lected such pair of values which is the best convenient for 
describing the observational cluster profile. Using the x 2 test 
we found the area of values po and r s with 68% confidence 
level. 



3.1.2 Stellar component 

Our sample includes clusters, which have a single, optically 
dominant cD galaxy near their centers, and the mass of the 
stars in the centra l galaxy was accounted through the Jaffe 
model (|jaffeHl983T ) with same parameters as in SA07. We 
assumed, as SA07 a mass of 1.14 x 10 12 M for the stellar 
mass of each central galaxy in the sample. 



3.1.1 Modeling of hot gas emission 

After 3D reconstruction of gas density field of pair, we were 
able to obtain the surface brightness in X-ray energy band. 
The emission flux, that comes from the observed regions, is 
proportional to the square of the integral density of hot gas, 
p 2 , along the line of sight. We used the metal abundance of 
0.3 and took n e /n p = 1.17 and p g = 1.35m p n p , where n e 
and n v are the electron and proton concentrations, respec- 
tively (|Vikhlinin et all \l99$ )). We have to note that from 



3.1.3 Best-fit 

All computations were made using x 2 -test 

x =>.( ) » ( 14 ) 

V a ° bs 

where Vobs and V mo de.i are the observed and modeled values, 
respectively. The regions with residual substructure, noticed 
by SA07, were down-weighted in the mass analysis, similarly 
to SA07. 
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3.1.4 Mass profiles scaling 

Density profiles of DM are approximatively universal on a 
large mass scale (Navarro, Frenk, & White 1995; Navarro, 
Frenk, & White 1996)0 

In our work we scaled mass profiles with -R200 an d M200 
(here -R200 is the radius within which the mean halo density 
is 200 times the critical density, and M200 is the total mass 
of galaxy cluster inside a sphere with radius i?2oo)- 

The total mass within a given overdensity A is defined 
in the present work as 



M A 



4 3 

— irApcriticRAi 



(15) 



where, as seen in the Introduction, R A = cat's is the radius 
within which the mean cluster overdensity is A times p cr itic 
and the relation with the concentration ca and the scale 
radius r a holds by definition of the NFW mass profile. In 
the present paper, we assumed A = 200. 

An important issue that we are dealing with in this 
paper, is to understand if the baryonic content in the clusters 
has influences on the slope of the density profile. To this 
aim, for each object we calculated the fraction Mb/M2oo- 
The quoted ratio Mt,/M2oo will be obtained subtracting from 
the total mass the DM mass in each cluster. So, the baryonic 
fraction will be derived as 



M„ = 



800n n \R 3 R 3 1 

Pcritic 1-TL200, DM + b ™ ^200,DM] 



(16) 



We will then plot the obtained values of a, namely the 
fit to DM component through a gNFW model, versus the 
baryonic content Mb- 



4 RESULTS 

In this section we show and discuss the results for the to- 
tal mass profiles, the DM profiles, the mass-concentration 
relation, and the inner slope of the DM profile. 



4.1 Total mass profile 

In Tab. [2] we summarized the results obtained from mod- 
eling the total mass profiles with a NFW model (a = 1), 
similarly to SA07. We plotted just the result concerning the 
NFW fit and discard the singular isothermal sphere (SI) 
(p(r) oc r~ 2 ) since the result of the comparison to the total 
mass profile with a NFW mode and that with a SI is simi- 
lar to SA07. The first column includes the names of galaxy 
cluster studied, the second column the parameter r s , the 
third C200 = R20o/r s , the fourth shows the value of -R200, 
the fifth includes the value of total mass at -R200, and the 
sixth shows the goodness of fit (\ 2 /d.o.i.). Confidence limits 
are 68 per cent. For most clusters the NFW model is a good 
fit, similarly to SA07. 




_J I I L_ 

10 100 
M eo[ 10 14 M Su „ 



Figure 3. The mass-concentration relation for the total mass 
(C200-M200). Solid line: best-fitting model with co = 6.83 ± 0.64 
and a = -0.19 ± 0.07 (6 = 1 fixed) for model [19] Dashed line: 
best-fitting model with c = 7.73 ± 1.07, a = -0.56 ± 0.15 and 
b = 0.80 ± 0.25 for model l20l (both have 68 per cent confidence 
level) . 



4.2 Dark matter profile 

In Tab. we summarize the result from the analysis in 
which cluster mass distributions were separated in DM halo 
(fitted with a gNFW model), the X-ray emitting halo, and 
the luminous mass of the BCG. 

We considered as SA07 two cases. In the first case, we 
examined models in which the inner slope of the DM profile 
was fixed at a — 1, corresponding to the NFW model, and 
in the second, we decomposed the total mass in its DM, 
diffused gas, and BCG, component (see Sect. 4.3). Again, 
except some cases, the NFW model provides a good fit to 
the DM profiles in clusters. 

We also calculated the concentration parameter for each 
cluster and this was used to obtain the mass-concentration 
parameter. This was done calculating cw = C200 and 
M vir = M2ocQ where 

M200 = ^TT200p crlU cRloo (18) 
As in SA07, we define the concentration parameters as 

r . — f -total I darfc-|8| 
Cvvr — {I v ir / 1 s JlS 

The mass-concentration relation measured from the 
Chandra, shown in Fig. [3] shows that the mass- 
concentration relation decreases when mass increases. The 
simplest analytic form that des cribe the mass-con centration 
relation is a power-law model |Dolag et al.ll20o3 ) 



c(z) 



M 



co 



1 + 2 V 8 x lO^ft-Wg 



(19) 



7 Note that in SA07 c v i r , and M„; r are connected through 



As previously reported, more r ecent simulations (e.g 



Power et all |2003|. lHavashi et al.|l2004 and lNavarro et al.|[20o3. 
Navarro et al. I I2OO4 IStadel et al. 2009) showed that density pro- 
files are better fitted by the Einasto profile. 



-TTr vir A c (z)p cr iti. 



(17) 



IShaw et al] ||2006|') . with A c = 178£2 m (z) - 45 llLahav et al.lll99lll 
8 Both virial radii and virial masses are calculated for the total 
mass model, including all mass components. 
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Table 2. The list of results from modeling the total mass profiles with Navarro- Frenk- White model with slope a = 1. First 
column includes names of galaxy cluster, second shows the r s parameter from NFW model, third shows the concentration 
parameter C200 = R200/RS for the NFW model, fourth shows the radius within which the matter density is 200 times the 
critical density for total and dark matter profiles (in Mpc), fifth and sixth demonstrate the total mass in clusters at radius 
-R200 and xVd-o.f. goodness of fit. Confidence limits are 68 per cent. 
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Table 2. Continued. 



Total mass (dark mattcr+gas+galaxics) 
NFW profile (a=l) 
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2 
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1.88 
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Table 2. Continued. 



Total mass (dark mattcr+gas+galaxics) 
NFW profile (e»=l) 

r s C200 R200 M200 x 2 /d-o.f. 

Mpc Mpc 10 14 M 



RXCJ2014. 8-2430 





50 


+0 
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07 
04 
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88 


+ 
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54 
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1 


95 


+ 
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2 


83 
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33 
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For example, I Shaw et al.l (|2006T ) by using the previous 
relation, with their simulated data, found Co = 6.47± 0.03 
and a = -0.12r jr. 0.03 for 68 pe r cen t confidence limits i n 
agreement with lEke et~al] i|200ll ). and lBullock et al.l l|200ll) . 
In a recent work by Kulinich et al. (2012), a new method to 
calculate a halo concentration parameter, and taking into 
account the halo overdensity and merging, was proposed. 
In our research, we used also a power-law model, obtaining 
co = 6.83± 0.64 and a = -0.19± 0.07 (6=1 fixed). As we 
show in Fig. [3] using the measurements obtained with our 
method above, the relation between concentration and total 
masses for CDM halos is represented not so well from Dolag 
et al. (2004) formula (Eq. 22). In agreement with SA07, 
Fig. 1, the model with 6=1, provides a poor fit for the 
observations, lying systematically below the data at lower 
masses and above it in the highest mass range. 

In order to improve the descripti on of Chandra data , 
we assumed, as SA07, differently from iBullock et all (|200ll ) 
to model a redshift evolution of the form (1 + z)~ b , with 6 
free 



C() 



M200 



(1 + z) b \8 x lO^ft-iMo 



(20) 



The motivation to the new redshift dependence comes from 
the fact that cl usters are more com plicated objects than 
those assumed in lBullock et al] (|200lh (see discussion in Sec. 
3.2 of Sa07, and Zhao et al. 2003). 

The results from a fit using Eq. (|20[) . with Co, a and 6 all 
free, are shown by the dashed line in Fig. [3] This model pro- 
vides an improved description of the data with \ 2 =1-3 and 
best-fitting parameteres Co =7.73± 1.07, a = -0.56 ± 0.15 
and 6 =0.80± 0.25 (68 per cent confidence limits). Note, that 
RXJ0848. 7+4456 and RXJ0849+4452 have been excluded 
from the fitting, since they give small values for concentra- 




M b 10"M Sun 

Figure 4. ci-Mj, relation for all sample of galaxy clusters. The 
data were obtained fitting the DM profiles of clusters with the 
gNFW model. The baryonic mass was obtained, as described in 
the text, by subtracting the DM component from the total mass. 



tion parameter (C200, 0.821q;qi and 1.141q io> respectively). 



4.3 The inner slope of DM profiles 

Finally, we fitted the DM profiles by means of a gNFW 
model to have hints on the inner slopes of the density pro- 
files. The results are summarized in Tab. [3] The value of a 
obtained is a = 0.94 ± 0.13 (68 per cent confidence limits). 
This result is consistent with CDM predictions. 

The values of the slope are in the range 0.5 < a < 1.8 
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and for the largest number of clusters the density profile is 
compatible with the expectation from ACDM model 0.7 < 
a < 1.5 (SA07). Only some clusters have slopes smaller than 
0.7 and they are usually characterized by an higher Mb with 
respect to clusters with larger a. 

4.4 Inner slope and baryons content 

As discussed in the introduction not all clusters have pro- 



files best fitted by the NFW model (e.g., ISand et al.l 


2002. 


Sand et al. 


12004 Sand et all 120081; iNewman et al l 


2009. 


Newman et al. 


2011). Among the processes able to pro- 



duce flat profile, the role of baryons has been severa l time s 
claimed (see the introduction). Recently. iDel Popolol (|2012f ) 
studied the effect of baryons, and the central BCG in shap- 
ing the inner density profile. It was also shown how higher 
content of baryons can give rise to profiles similar to A611 
and A383 which have slopes flatter than NFW profile, and 
small baryonic content can produce steeper inner profiles 
in agreement with the NFW model (MACS J1423.8+2404, 
RXJ1133). 

Motivated by the previous arguments, we studied the 
eventual correlation among the inner density profile slope, 
a, and the baryonic content of the clusters that we studied. 

Similarly to SA07, we performed an analysis in which 
the total mass was separated in the DM halo, the X-ray 
emitting halo, and the luminous mass of the BCG (this is 
the second case mentioned in Sect. 4.2). We fitted the DM 
profile with a gNFW model and calculated a, and reported 
the values in Tab. [3] From the analysis the baryonic content, 
Mb is calculated summing the diffuse gas content and the 
baryons contained in the central cD galaxy in each cluster. 
In Fig. |31 we plot the inner slope, a, versus the baryonic 
content at the radius of -R200 ■ A2255 has been excluded dur- 
ing fitting, because this galaxy cluster gave small uncertain 
value of x 2 - The plot shows a tight correlation between the 
inner slope a and the baryonic content: the larger is the 
baryonic content the flatter is the inner profile. This results 
is in agreement with studies of the role of baryons on the 
inner slope of clusters (e.g.. IDel Popolol 120091. IDel Popolol 



l2012l , lRomano-Dfaz et~aUl2008l , iRomano-Dfaz et aUfeOOgT 
An important issue to stress is the fact that we assumed, 
as SA07 that all the central BCG has the same mass, so that 
the difference of the baryon content from cluster to cluster is 
connected to the diffused gas content. The expectation that 
the central baryonic content is the only or the fundamental 
contribution to shape the inner density profile is only partly 
true. To understand this, it is interesting to discuss more 
in detail the role of the total baryonic mass, Mb, and the 
central one. 

Baryons and the baryon diffused component pres- 
ence produce s in general a fla t tening of the den- 
sity profile dEl-Zant et alj l200ll. lEl-Zant et al.1 |2004 
Romano-Diaz et al.ll2008l,|Pel Popololl2009l.lGovernato et all 



2010T ). and rounder halos than those seen in DM simu- 



lations (jDebattista et al.ll2008l ; lAbadi et al.l l201dh . The fi- 
nal configuration of a cluster is fixed by the initial quan- 
tity of baryons present in the proto cluster and by col- 
lapse/formation process. Then it is logic to expect that 
the final central baryonic content and the BCG mass 
is somehow correlated with baryonic and total cluster 
mass, M c i. For example, IWhilev et al.l (|2008l ) . found that 



Mbcg oc M®i 4 or M^' 5 depending on the feedback model 
used. Also a correlation between BCG luminosity and 
cluster X-ray luminosity was found by several authors 
dSchombertl (l988l; lEdge fe Stewart! Il99ll; lEdge fc Steward 



Il99ll ; iHudson fe Ebelind 1 1997ft . IWhilev et al.l d2008 ) 



sured the quoted correlation as Mbcg oc M 



mea 

0.12±0.03 



for 



K band ma gnitudes inside a dia meter of 37 kpc (radius of 



id mag nitudes ins ide a dia meter ot 
kpc). lBrough et al.l (|2008h found L 



oc M 



0.11±0.10 



13/i Kpo. iorougn et aij i|zuua i rouna Ijbcg <-»- j^ c ; 
at K band inside V2h~ x kpc (several other results are given 



Lin fe Mohr|l2004l; iPopesso et all 120071 1 Yang et al.H2008l ; 
1] 



lHaarsma etaU |2010| ; Fedeli 2012). The flattening of the 
slope can be interpreted as due to the fact that the pres- 
ence of a larger quantity of baryons guarantees a larger 
transfer of energy and angular momentum from baryons to 
DM, with the result th at DM moves to larger orbits, reduc- 
ing the inner dens i ty (IDel Popolol 120091. |Pel Popolol |20l3: 
El-Zant et~aT]|200ll.lEl-Zant et al.ll2004l ; iRomano-Pi'az et al l 
20081 . iRomano-Pi'az et al.ll2009l ). 



As already reported, NFW N-body simulations pre- 
dicted a profile with inner slope oc r^ 1 and an outer one 
oc r -3 . More recently, it has been shown that the density 
profiles are better fitted by an Einasto profile which becomes 
shallo wer towards the centre of the halo (e.g. iNavarro et al.l 
[2010J). In any case, the lower values for the inner slope 
obtained in N-bod y simulations are a = 0.8 at 120 pc 
l|Stadel et al.l l2009) in agreement to ana lytical and numeri- 
cal works to solve the Jeans equation (|Austin et al 1 l2005l ; 
iPehnen fc McLaughlin! [2005) . Inner density slopes in the 
range 0.7 < a < 1.5 could be considered in agreement with 
ACPM predictions (SA07). 

However, as reported in the introduction, slopes flatter 
than those obtained in the simulations, a < 0.7 are ob- 
served in some clusters llSand et al. 2002, ISand et al. | |2004 
ISand et al.ll2008l ; INewman et al.ll2009l . INewman et~ai1l201ll ) 
and the largest numbers of dwarf galaxies and LSBs, and, 
as previously reported, also some of our clusters have slopes 
a < 0.7. 



This discrepancy does not nece ssarily imply a problem 
for the ACPM model (e.g. see PP09; iGovernato et alfeoiOT ) . 



Given the several and noteworthy pieces of evidence 
supporting ACPM on large scales, the discrepancy could be 
connected to the fact that baryonic physics is a fundamen- 
tal issue in cluster formation. It could originate just because 
we are comparing dissipationless systems generated by N- 
body simulations with real, dissipational structures, whose 
physics is different from the dissipationless physics typical of 
PM. We should not forget tha t the inner 10 kpc of clusters 
are dominated by baryons (e.g., Sand et aLll20 02, San d et al 
20041. ISand et al.ll2008l ; INewman et al.ll2009l , INewman et al 



20111 ) whose presence strongly influence the PM distribu- 
tion. While baryons can steepen the inner slope of the 
density profil e of clusters through th e adiabatic contrac- 
tion of PM (iBlumenthal et all Il98rj ; iGnedin et al] |2004| ; 
iGustafsson et alT 20061 ). heating of PM due to dynamical 
friction with cluster galaxies can counteract the adiabatic 
contraction effect and flatten the inner profile ( El-Zant et alj 
l200ll , lEl-Zant et aT]|2004l ; IRomano-Pi'az et al.il200Sh . 
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Table 3. Results on dark matter profiles for the clusters, from the fits that includes DM, X-ray emitting gas, and BCG. 
The first column includes names of galaxy cluster, while the second, third and sixth column refers to the NFW fit: second 
shows the r s parameter from NFW model for dark matter, third shows the concentration parameter C200 = ^200 /Rs , sixt the 
X 2 /d.o.f. goodness of fit. In the case of the gNFW, we give the inner value of a, in seventh column, while in eighth column 
is reported the xVd.o.f. goodness of fit. Fourth column shows the mass of baryonic matter in clusters, and fifth shows the 
fraction of baryonic matter. Confidence limits are 68 per cent. 



Dark matter profiles 

NFW profile (a — 1) NFW (a is free parameter) 
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4 21+ - 55 

^■ z -0.88 


20.46l 3 ;^ 


0. 


105 


1 


,05 


1 


,041° 


45 
28 


2 


.13 


A1835 


0. 


■87l° 


18 
27 


2.36l°' 3 4 


41 93+ 5 - 16 


0. 


183 


1 


.83 





,68l° 


13 
.15 


1 


.77 


A1914 


3 


.051° 


48 
37 


74+0-36 


93.76l° 6 27 


0. 


121 


1 


,21 





■43l°„ 


18 
,18 


1 


.28 


A1995 





88l° 


33 
28 


9 77 +0.38 
' ' -0.33 


56.22l 6 4 -I 


0. 


124 


1 


24 





,551° 


18 
.18 


1 


.22 


A2029 





.37l° 


17 

.19 


6 39+0- 41 
o.oy_ 43 


42 42+ 5 ' 18 


0. 


126 


1 


,26 





.911° 


37 
,17 


2 


.20 


A2034 





■90l° - 


23 
.22 


9 , 7 +0.49 
' -0.42 


91 99 + 1-93 

zi.zz_ 2 18 


0. 


106 


1 


.06 


1 


,06l° 


37 
,27 


2 


.38 


A2052 





13+°' 


03 
04 


9.63±\i 9 3 


13.45l 2 ; 46 8 


0. 


147 


1 


,47 


1 


,62l° 


27 
.22 


2 


.32 


A2063 





22+°' 


11 
05 


b - 89 -0.43 


20.061 2 ' 47 


0. 


194 


1 


,94 





,80l° 


27 

.28 





.92 


A2065 


1 


99 + 0.21 
zz 0.18 


2 02+ - 34 

z - uz -0.46 


94.05H?; 27 


0. 


238 


2 


,38 





41+0 

^^-0 


29 
.27 


1 


.30 


A2124 


0. 




04 
02 


1 n 7 ^+l-84 

10.761! .gi 


20.67i 2 ;i 5 9 


0. 


157 


1 


,57 





-9ll°o 


12 
16 


2 


.47 


A2142 


0. 


.501° 


12 
11 


5 in+°- 43 


2l.69l 3 ; 2 8 8 


0. 


072 





,72 





,931° 


12 
17 


1 


.29 


A2147 





.161° 


02 
03 


10.061} 84 


18.09+ 2 ' 27 


0. 


165 


1 


.65 


1 


-12l° 


19 
19 


2 


.28 


A2163 


2 


ool° 


47 
45 


1 ss+°- 27 

1.8S_ 23 




0. 


129 


1 


.29 





■Bit?' 


14 
.15 





.75 


A2199 


0. 


■24l° - 


09 
10 


r qo+0.45 


24.25l 2 ;?| 


0. 


195 


1 


,95 





90l°o 


39 
.10 


2 


.17 


A2204 


0. 


.761° 


21 
17 


2 56+°' 42 

Z.OD_ 54 


43.67l 5 -; 82 9 


0. 


186 


1 


,86 





,85t° 


36 
37 


2 


.28 


A2218 


0. 


.211° 


05 
03 


D - 1,5 -0.65 


8-95l 2 ; 9 5 3 


0. 


138 


1 


,38 


1 


,78t° 


14 
.15 


2 


.10 


A2219 


0. 


.661° 


15 
17 


3 24+°' 43 

-0.53 


47 06+ 5 ' 38 

11 .uu_ 4 19 


0. 


128 


1 


,28 





,83l° 


37 

.33 


1 


.18 


A2244 


1. 


.121° 


25 
.27 


9 90 + O.34 
z - zo -0.31 


35.20_3 g 2 


0. 


096 





,96 





-BUS 


16 
.17 


2 


.19 


A2256 


2 


67l° 


44 
17 


1 06+ ' 47 

r.uu_ 37 


H7 77 +5.97 
' • ' ' -6.92 


0. 


120 


1 


,20 


1 


-12tS 


11 
27 


1 


.11 
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Table 3. Continued. 



Dark matter profiles 

NFW profile (ct=l) NFW (a is free parameter) 

Name r s c 2 oo M b fM b /M tot X 2 /d.o.f. a xVd.o.f. 
Mpc 1O 13 M 



A2319 


2 49+ 22 


1 08+°' 37 

r.uo_ 32 


71 45+ 8 - 14 


0.155 


1 


.55 


o.sei 9 :! 2 


2.22 


A2390 


, -, 7 +0.34 
J " - 1 ' -0.26 


1 91 + ° 42 
-"-•^-0.32 


50 S9+ 5 67 
du - oy -4 .86 


0.180 


1 


.80 


0.94l»:l 2 


2.02 


A2462 


u - lo -0.06 


12.63l°o;f 7 


28.331 3 ;! 4 


0.136 


1 


36 


i 91 +0.15 


2.05 


A2537 


40+ 012 


4.i8l° : 3 4 4 4 


31.45l 4 ;f 9 


0.195 


1 


.95 


1 18+ ' 13 

±i±o -0.11 


2.66 


A2589 


n ,,+0.09 
u - ol -0.07 


5.27l :S 


18.43l 3 ;^ 


0.183 


1 


.83 


i.34i° :l 7 


2.48 


A2597 


n ,o+o.05 

U - i8 -0.03 




12.3411; 27 


0.114 


1 


.14 


0.92l° :l 8 


1.12 


A2634 


15+ " 5 




99 i i +2 .63 

zz.ii_ 2 18 


0.186 


1 


.86 


-1 9 ,+0.19 
1 - zl -0.24 


2.47 


A2657 


33+° 11 


5 14+ 043 


19.531 2 -\% 


0.197 


1 


.97 


-1 -IO + 0.13 


2.55 


A2667 


1 00+ 023 

±uu -0.12 


„ „, +0.47 
z ' ZJ --0.25 


37.591 4 ;! 2 


0.190 


1 


.90 


74+ - 14 

u - '^-0.11 


2.28 


A2670 


41+ 013 


5 m 4 " 51 

°- ui _0.59 


30.09l 3 ^ 


0.180 


1 


.80 


74+ 12 

u - '^-0.12 


1.21 


A2717 


UJU _0.04 


5.88121 4 


25.131 2 ; 34 


0.194 


1 


.94 


i n3+ 014 


2.83 


A2744 


i ooO 30 
1 - OZ -0.37 


1 52+ - 35 

1 - OZ -0.57 


™ ,! " -4.72 


0.220 


2 


.20 


o.7oi° :l 2 3 


2.36 


A3112 


21+ 04 

UiZ1 -0.03 


9 16+° 43 

^'^ -0.63 


13.6ll 2 - 27 


0.072 





.72 


0.83l° :ll 


1.27 


A3158 


1 20+ 013 


9 n9 + 0.60 
z ' uz -0.61 


16.76l 2 - 27 


0.065 





.65 


o.95i :i 4 


1.22 


A3266 


n ,,,+0.05 


„ 1n + 0.53 


42.05l 4 ;f 7 


0.126 


1 


.26 


1.2312:1? 


1.21 


A3376 


n , 9 + 0.03 
UJZ -0.04 


7 01 + li11 


30.671 4 ; 27 


0.115 


1 


.15 


, fi o + 0.35 


2.29 


A3391 


21+ 05 

uzl -0.03 


9 06+ 54 

a ' UD -0.76 


30.6ll 4 ' 2 27 


0.161 


1 


.61 


o.72i° :l 9 


2.02 


A3395 


71+ 011 

u - ' 1 -0.20 


3 34+ 45 

°-° -0.47 


45.40l 6 5 s 2 7 


0.134 


1 


.34 


0.82l :l 4 


2.07 


A3526 


10+° 04 

U.1U_ 03 


lO.lSll;^ 


18.241 2 %% 


0.220 


2 


.20 


0.69l :l 4 


2.66 


A3558 


22+ 04 


7.89l 9 o;f 3 


28.64l 4 ; 2 7 


0.134 


1 


.34 


o.76i° :l 8 8 


1.33 


A3562 


22+ ' 05 

U - ZZ -0.06 


7.96l 9 ;f 7 


33.07+ 4 ; 27 


0.235 


2 


.35 


67+° 18 

UiUI -0.19 


1.17 


A3571 


25+ 06 

u - zo _0 03 


8.52±°-» 


39.35_ 3 2S 


0.129 


1 


.29 


0-9lloT9 


2.36 


A3667 


n 9 , + 0.06 
uzo -0.04 


8-86l° ; 6 6o 


44 18+4.56 


0.164 


1 


.64 


o.87i° :l 7 


2.22 


A3827 


U - Si -0.23 


4 l, + 47 
^• lo -0.57 


70 7Q + 8 - 36 


0.089 





.89 


n ro + 0.24 
UOO -0.34 


2.02 


A3921 


n 46+ - 08 


4.87l°o; 4 31 


18.56l 2 ; 3 74 


0.062 





.62 


1 16+ - 48 

1 - 1D -0.29 


2.33 


A4038 


15+ ' 02 
u - 10 -o.oi 


8.78ll;° 2 3 


14.98ll; 4 ° 8 


0.180 


1 


.80 


0.77l° :l 2 3 


2.32 


A4059 


17+ - 02 

u1 ' -0.00 


8.86ll;°o? 


13.64ll;H 


0.127 


1 


.27 


0.83l° :l 2 5 


2.37 


AWM4 


15 +0 03 
u - lo -0.02 


6.61 + 1' 32 


IK ci +1.47 

15.61IJ 47 


0.238 


2 


.38 


, - 9 +0.25 
1 - oz -0.22 


2.22 


CLJ1226. 9+3332 


1 54+ 013 
o.n 


1 54+ ' 42 

1 '°*-0.34 


94.12li° 27 6 


0.065 





.65 


0.5212:;? 


1.45 


Coma (A1656) 


2.27lg; 33 


1 26 + ° 35 


76.82l 7 ; 6 5 2 


0.183 


1 


.83 


o.63i2:l 6 6 


1.24 


IIIZw54 


12+ 001 

u - lz -0.02 


10-6511 2 | 


17.99l 2 - 74 


0.204 


2 


.04 


i.43i2: 2 2 ? 


1.16 


ISCS J1438. 1+3414 


6 41+ 35 

°-* 1 -0.42 


43+ 05 


316.22ll 4 2 8 92 


0.107 


1 


.07 


o.35l2:l 4 4 


2.73 


MACSJ0011. 7-1523 


n cc+0-23 
°- 65 -0.13 


, tjo + 0.33 


59 Q7+ 4 - 88 


0.124 


1 


.24 


o.esig-.ii 


1.33 


MACSJ0159. 8-0849 


44+ 011 

0.11 


4 85+ 022 

^• Od -0.32 


52.57l 4 :f 2 


0.119 


1 


.19 


g 79 +o.24 

u - ' z -0.36 


1.19 


MACSJ0242. 6-2132 


n 97 +0.06 
u - z, -0.04 


7 ,o + 1.24 
'■ 18 -1.13 


48.48ll: 2 « 


0.176 


1 


.76 


64+° 17 
u.u<±_ 17 


2.33 


MACSJ0429. 6-0253 


55+ ' 23 


2 Q6+ 1 ' 36 
zaD -1.15 


31.31I 3 : 24 


0.169 


1 


.69 


93+° 17 


2.38 


MACSJ0647. 7+7015 


, ,„ + 0.15 
°- 1 -0.16 


90+ 011 

U.»U_ 21 


98.25l?- 2 8 


0.109 


1 


.09 


77+ ' 37 

u - ' ' -0.27 


2.12 


MACSJ0744. 8+3927 


3 14+ 019 


n oi +0.27 

0.81_ 16 


97-8ll 9 7-.99 


0.091 





.91 


0.65_ 19 


2.28 


MACSJ1115. 8+0129 


1 00+ 019 

luu -0.12 


, no + 0.33 
l- 8d -0.25 


49.381 4 :™ 


0.209 


2 


.09 


64+ - 26 


2.21 


MACSJ1311. 0-0311 


50+ - 07 

u - ou -0.06 


4 51+ 23 

^• ol -0.15 


68.15l 7 :l 5 2 


0.123 


1 


.23 


o.54i2:f 3 


1.18 


MACSJ1423. 8+2404 




rj Q q + . 3 5 

z ' yo -0.23 


34.90l 32 5 


0.072 





.72 


o.63i2:i 8 8 


2.88 


MACSJ1427. 6-2521 


26+ 07 
u - ZD -0.02 


7.58l 9 ; 77 


34.75l 3 - 3 7 


0.138 


1 


.38 


o.84i2: 3 2 3 


1.18 


MACSJ1720. 3+3536 




4.85l°o: 3 46 


37.6ll 4 4 : 29 


0.092 





.92 


79+° 17 

u - ' "-0.17 


2.72 


MACSJ1931. 8-2635 


0.6712:15 


3.61±°;£ 


42.24l 71 18 


0.129 


1 


.29 


n 77 +0.38 
u - ' ' -0.33 


1.44 


MACSJ2129. 4-0741 


1 96+ - 22 

laD -0.37 


, o 9 + 0.28 


95.79l 7 6 .S 


0.099 





99 


o.64i2:l 8 8 


2.46 


MACSJ2229. 8-2756 


n 97 +0.06 
u - z ' -0.03 


7 64 +0 <37 

'• D ^-1.16 


35.641 4 82 


0.077 





.77 


o.84i2:; 9 9 


2.33 


MKW3s 


n ,0 + 0.03 
ulo -0.03 


I0.47ll; 7 l 


20.87+ 3 99 


0.173 


1 


.73 


o.85i2:l 7 


2.22 


MKW4 


o ri5+ 001 
u.uo_ 00 


18.93+llf 


7 69+ 1 94 

1 .u»_2.19 


0.346 


3 


.46 


7 +0.28 
1 -'°-0.27 


2.37 


MKW8 


10+ ' 02 


i4.68i;;i 4 


8 34+ 2 * 37 
- J ^-1.87 


0.111 


1 


.11 


, no + 0.18 

1.83_ 17 


0.92 


PKS0745-191 


n 00 + 0.03 
U - J8 -0.02 


- 7 c + 0.53 


30.04+ 2 1% 


0.089 





.89 


o.88i2:l 8 9 


0.82 


RXCJ0043. 4-2037 


14+ 003 


7 42+ 1 89 

' '^ z -1.03 


10.29l 2 3 ^ 


0.161 


1 


.61 


i-43l2'.; 8 8 


2.11 


RXCJ0232. 2-4420 


, 97 +0.12 
-0.13 


1 oS + 047 


16.5lll:f 5 


0.090 





.90 


, ,,+0.34 
1 - ol -0.3S 


2.22 


RXCJ0307.0-2840 


55+ 017 
u.oo_ 16 


9 - 9 + 0.38 
^•■^-0.37 


21.28+2 ll 


0.128 


1 


.68 


, 9 ,+0.19 


2.38 


RXCJ0516. 7-5430 


74+ 019 

Ui '^-0.16 


2 3^ + 55 




0.126 


1 


.26 


1 34+ 017 

1 -°*-0.16 


2.72 


RXCJ0547.6-3152 


41+ 06 


3 64+ 77 


18.82l|:l| 


0.150 


1 


.70 


1 64+ 015 


2.34 


RXCJ0605.8-3518 


34+ 02 

u - J ^-0.02 


3 74+ - 23 

°' ' -0.33 


13.8411H 


0.117 


2 


.17 


1 54+ ' 19 

± - d *-0.19 


2.38 


RXCJ1131. 9-1955 


U - '°-0.13 


2 25+ 114 


20.92+ 2 28 


0.115 


1 


.35 


1 42+ 018 


2.88 


RXCJ2014. 8-2430 


46+ 07 

u -™-0.04 


o co + 0.34 
3.t)8_ 47 


18.27llj s 3 


0.123 


1 


.73 


1 64+ 019 


1.83 


RXCJ2129.6+0005 


42+ 07 


3 17+ 9 ' 48 
°- L ' -0.34 


1 2 97 +!.33 
iz ' z '-1.26 


0.115 


1 


.45 


, 7 „ + 0-24 
1 - '°-0.22 


1.33 
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Table 3. Continued. 



Dark matter profiles 

NFW profile (a — 1) NFW (a is free parameter) 

r s c 2 oo M b fM b /M tot x 2 /d.o.f. a xVd.o.f. 

Mpc 1O 13 M 



RXCJ2337.6+0016 


29+ - 06 

u -^ a -0.03 


4 y 2 + 0,55 


13.T5+| : || 


0.131 


1.61 


1 67+ - 19 

± ul -0.19 


1.18 


ZwCL1215 


41+ 08 

"■^-0.16 


E i 9 +0.42 


36.33_ 2 3g 


0.118 


1.08 


74+0- 14 

u - '^-0.13 


1.14 


RCSJ0224-0002 


51+ 009 
U0i -0.17 


9 AS + 35 
2 - 48 -0.35 


22 45+0- 19 


0.113 


2.13 


1 02+ 015 

1 - uz -0.14 


1.39 


RCSJ0439-2904 


20+ 01 


, Qrj + 0.35 




0.207 


2.37 


1 84+ 036 


2.20 


RCSJ 1107-0523 


34+ 03 


n 7C+0 36 
z - ,o -0.35 


1 9 KK + 1-47 


0.219 


2.39 


, -O+0.54 
1 - oo -0.53 


1.25 


RCSJ1419. 2+5326 


n „„+0.03 
u - zz -0.05 


5.74t° ; 5 3 3 


33.91 + 3 * 7 


0.222 


2.62 


67+ 012 

u -°' -0.13 


1.28 


RCSJ1620+2929 


31+° 15 


3.68±S;S 


14.16±ii£ 


0.147 


1.87 




1.17 


RCS2156+0123 


, --+0.22 
L - 1 ' -0.34 




42.97±|;g 


0.181 


1.88 


0.58ig£ 


1.25 


RCSJ2318+0034 




, n „ + 0.29 

1.08_ 27 


51.45+. 5 - 


0.141 


1.81 


64+ 015 

u -"^-0.14 


1.10 


RCSJ2319+0038 


U -°'-0.09 


9 Q9 + 0-36 


Q1 Q1+ 8 - 88 


0.204 


2.24 


54+ 015 

U - J *-0.16 


2.36 


RXJ0439. 0+0520 


i 7 +o.05 

u ± ' -0.03 


7 o-i +0.27 
'■•5 1 _0.77 


11 43+ 216 

.37 


0.225 


1.25 


1 04+ 029 


2.46 


RXJ0848. 7+4456 


, OO+0.34 

i.oo_ 27 


57+ - 15 




0.241 


2.81 


1 01+ 014 

'■■"'--0.15 


1.66 


RXJ0849+4452 


, 7 , + 0.37 
1 - '°-0.25 


74+°' 21 

u - '*-0.12 


67.87± 6 6 .f 4 


0.165 


1.45 


uoo -0.16 


1.28 


RXJ0910+5422 


S3+ 015 


2 44+0- 15 


66.84+.?^ 


0.143 


1.03 


o.63+° .f 3 


2.46 


RXJ1113. 1-2615 


26+ 06 
U - /D _0.04 


2.98+. ,' 23 


12 72+ 1 - 49 


0.219 


2.39 


i-04l° ; 33 7 


2.22 


RXJ1221. 4+4918 


1 01+° 15 


9 n , + 0.27 


75.08+. S; 32 


0.151 


1.61 




1.27 



5 DISCUSSION 

As discussed in the paper, our analysis used Chan- 
dra data and similar methods to that of SA07, but 
our sample is larger. Confronting our tables and those 
of SA07, we see several clusters common to the two 
studies (e.g., A1795; A2029; A478; A1413; A2204; 
A383; A963; A1835; A611; A2537; MACSJ0242.6-2132; 
MACSJ1427.6-2521; MACSJ2229.8-2756; MACS J 1391.8- 
2635; MACSJ1115.8+0129; MACSJ1720.3+3536; 
MACSJ1311. 0-0311; MACSJ1423.8+2404). The inner 
slopes of these clusters are compatible in our and SA07 
study. 

Other clusters, like A383, A611, A963, A1835, A2029, 
A2204, A2589 were st udied by others autho rs. For example, 
A611 was studied by iNewman et all (|2009h who combined 
weak lensing from multicolor Subaru imaging, strong lens- 
ing (Hubble Space Telescope) and stellar velocity dispersion 
measures (Keck Telescope), sampling the dark matter pro- 
file from 3 kpc to 3.25 Mpc. Newman et al. (2009) found 



values of r B 



320 



+240 



kpc and c = 5.lli'g (in agreement 



with SA07), but a < 0.3(< 0.56, < 0.65) at 68 per cent (95 
per cent, 99 per cent) CL. SA07 found a value of the slope 
A383 was studied by S08 finding a flat inner 
= 0.45l^2 5 ), and by Newman et al. (2011) 



a = 0.64 +u - M 
DM slope (a 

who found a < 1 at (95 per cent confidence) and a best 
fit (inferred from weak and strong lensing, kinematics and 
X-ray data) of a = 0.59±g"|§. SA07 found a < 0.8, we found 
a = 1.58t°;g. 

A fla t slop e wa s obtained also for A 9 63 b y 
Sand et all (120081), bu t iBartelmann fc Meneghettl (|2004h . 
Dalai fc Keetonl (|2003T ) and SA07 found values consistent 
with the NFW model . We f ound a = 1.52" 3 12 . In the case 
of A2029. lLewis et al.1 (|2003l ') found a = 1.19±0.04 similarly 
to SA0 7, and we got 0.9 1 +^. In the case of A1835, A2029, 
A2204, lArabadiis et al.1 \2004) by means of Chandra data 



obtained a ~ 0.9, a ~ 1.85, and a ~ 1.8, respectively, much 
larger than SA07 v alues and ours. 

As reported, INewman" et al.1 l|2009D . H ewman et al.l 
|201ll ). presented a detailed analysis of DM and baryonic 
distributions in A611, and A383 combining weak lensing, 
strong lensing and stellar velocity dispersion for the BCG, 
and findi ng slopes flatter than the NFW predictions . Ac- 
cording to lNewman et al.1 |2009l '). INewman et"ai1 l|201lf ). de- 
generacies in constraining the DM profile can be broken only 
simultaneously using the three techniques. In reality, the X- 
ray observations alone give information on clusters on cluster 
structure in the range 500-50 kpc. At smaller radii, temper- 
ature determination is limited by instrumental resolution, 
substructure (SA07), and they are also limited by "cooling 
flows" presence and t he breaking of assump tion of hydro- 
static equilibriu m (see lArabadiis et al.ll2004h . 

So, if the INewman et al.1 (|2009h and INewman et al.l 

(2011) point of view is correct, for clusters containing larger 
quantity of baryons the inner profile may be flatter than X- 
ray observations. SA07 tried to understand how important 
was the ro l e of c e ntral BCG i n MS2 137.3-2353 (for which 
ISand et al.1 l|2002f l. ISand et al.1 l|2004h found a flat profile). 
Varying the ra t io M /Ly they observed, i n agr eement with 
INewman et al.1 (|2009l ) and lNewman et al.l (|201lh a flattening 
of the profile, but in any case the s lope was still compat i- 
ble with a NFW fit ( differently from INewman et al.1 l|2009l) . 
INewman et~al1 i|201ll) '). 

Other systematic uncertainties that could change the 
results is the presence of non-thermal pressure support (due 
to gas motions, cosmic rays or magnetic fields) , however for 
relaxed clusters, bulk and/or turbulent motions, if present, 
could cause ch anges in the mass measurements of 10-20 per 
cent accuracy (|Nagai et al.ll2007l ; iRasia et al.ll2006l). Similar 
effects are expected by magnetic pressure iDolag fc Schindlerl 
|2000l ). 

Another important issue to recall is the fact that in 
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the present paper, as in that of SA07 the X-ray data only 
extend to about r2soo (approximately a quarter to a third of 
the virial radius in most clusters) , so, for example, the claim 
that the mass-concentration relation is appropriate for the 
virial radii, is not totally true. 



6 CONCLUSIONS 

We presented the reconstruction of the total mass (dark 
matter, gas and luminous matter) from the Chandra ob- 
servations of 129 massive X-ray luminous galaxy clusters in 
the redshift range 0.01 - 1.4. We estimated the total (Mtot) 
mass within -R200 in the range 3 - 293 xlO 14 M . 

In order to estimate the fraction of dark matter and gas 
we have used the conditions of hydrostatic equilibrium and 
spherical symmetry. Similarly to SA07, we performed two 
different analysis: the first concerning the total mass of the 
clusters, and the second decomposing the total mass in its 
diffuse gas, DM, and BCG components. 

Similarly to SA07, the NFW gives a good fit to the total 
mass and DM distribution. We obtained a best-fitting result 
for the inner slope of the DM profile in the clusters a — 
0.94 ±0.13. We also obtained a mass-concentration relation 
c oc M a /(1 + z) b , with a = -0.56 ± 0.15 and b =0.80± 0.25 
(68 per cent confidence limits) in agreement with previous 
results and simulations. 

Finally, we showed that there is a tight correlation 
among the inner slope a and the baryonic mass, M b , content. 
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